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'^ : ABSTRACT 

^ . Clusters of galaxies at redshifts nearing one are of special importance since they 

C^^ , may be caught at the epoch of formation. At these high redshifts there are very few 

^vg , known clusters. We present follow-up ASCA, ROSAT HRI and Keck LRIS observations 

CN I of the cluster RXJ 1716.6+6708 which was discovered during the optical identification 

Qs ' of X-ray sources in the North Ecliptic Pole region of the ROSAT All-Sky Survey. At 

0^ ■ z=0.809, RXJ 1716.6+6708 is the second most distant X-ray selected cluster so far 

1-^ ■ published and the only one with a large number of spectroscopically determined cluster 

^^- member velocities. The optical morphology of RXJ 1716.6+6708 resembles an inverted 



S-shape filament with the X-rays coming from the midpoint of the filament. The 



c/3 , X-ray contours have an elongated shape that roughly coincides with the weak lensing 

' contours. The cluster has a low temperature, kT = 5.66'^Q^g keV, and a very high 

velocity dispersion aios =15221]^5q km s~^. While the temperature is commensurate 
X ' with its X-ray luminosity of (8.19+0.43) x 10"^^ h^o^ erg s"^ (2-10 keV rest frame), 

C^ ■ its velocity dispersion is much higher than expected from the a — Tx relationship of 

present-day clusters with comparable X-ray luminosity. RXJ 1716.6+6708 could be an 
example of a protocluster, where matter is flowing along filaments and the X-ray flux 
is maximum at the impact point of the colliding streams of matter. 
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1. Introduction 



In hierarchical theories of structure formation, clusters of galaxies form from the high peaks 
in the original density field, thus they provide crucial constraints on the shape, amplitude, and 
temporal evolution of the primordial mass fluctuation spectrum. Despite their importance, the 
statistics for the abundance of high-z (>0.5) clusters are poor since they are so difficult to locate. 
Optical surveys of distant galaxy clusters are well known to have serious statistical shortcomings 



such as effects of superpositions of unvirialized systems (see among others Frenk et al. 199C ; Van 



Haarlem et al. 1997). In addition, at high redshifts it becomes difficult to detect enhancements in 



the galaxy surface density against the overwhelming field galaxy population (see Oke et al. 1998 
for the success rate of optical high-z cluster searches) . 

One of the cleanest ways to avoid sample contamination is the selection of high-z clusters 
by means of their X-ray emission. X-ray surveys are sensitive enough to detect distant clusters. 



Examples include MS1137-I-66, at z=0.78 and MS1054— 03 at z=0.83 in the Medium Survey (Gioia 
|et al. 1990| and Gioia and Luppino 1994; see also the detailed study on MS1054— 03 by ponahue 
et al. 199 j ). Other X-ray surveys being conducted with ROSAT archive data are also finding 
distant clusters, for example the WARPS cluster at 0.83 (RXJ0152. 7-1357, Ebehng et al. 1998) 
and the clusters discovered in the RDCS survey ( |Rosati et al. 199^ ; Rosati et al. 199^ ). The 
RDCS holds the record for the most distant X-ray selected cluster known at z~1.26 (Rosati et 
al. 1999). These X-ray surveys are demonstrating that they can go as deep as optical searches in 
locating distant clusters. 

Other techniques are successfully isolating clusters. Clustering around radio galaxies at 
redshift equal or greater than 1 has been reported (see among others 3C 184, z=0.996, Deltorn 
et al. 19971; MRC0316-257, z=3.14, [LeFevre et al. 1996|; 3C 294, z=1.786, [Crawford aild| 



[Fabian 1996 ; 3C 324, z=1.206. Small and Dickinson 1995). A new promising technique which 
is currently providing interesting results is the search for distant clusters in the near-IR (GIG 
J0848-I-4453, z=1.27 from 8 cluster members, Stanford et al. 1997[ ). However, with the exception 
of the survey of nine high redshift optically-selected clusters by Oke and collaborators ( pke et al 



1998 , Lubin et al. 1998 , Postman et al. 1998 ), most of these distant structures have spectroscopic 



determinations for less than one dozen cluster members, at best. 

Our group at the Institute for Astronomy in Hawai'i has been involved for several years 
in the optical identification of all the sources found in 84.7 deg^ centered at the North Ecliptic 
Pole (NEP) of the ROSAT All-Sky Survey (RASS, [Voges et al. 1996| ). The NEP region is the 
deepest area of the RASS. We reported in Henry et al. (1997) the detection of a distant cluster, 
RXJ 1716.6-1-6708, at z=0.8 with only 33 net photons. This weak detection revealed a very 
interesting object. The spectroscopy in hand at the time of the Henry et al. (1997) paper was 
limited to only 12 cluster members but there was already an indication (see their Fig. 1) that 
the galaxies in the inverted S-shaped filament were all part of the X-ray cluster. We report here 
follow-up observations of RXJ 1716.6-1-6708 at optical and X-ray wavebands. 



Very high-redshift X-ray selected clusters, such as the one presented here, are important to 



study. Optically selected clusters tend to be under- luminous X-ray sources ( Bower et al. 1997 and 
reference therein), an indication that generally they are not massive objects. As more accurate 
observing capabilities become available, the number of distant known clusters will steadily 
increase, and the counts may severely constrain cosmological models. However, we are still far 
from having large complete samples of high-z clusters to test the predictions of cosmic structure 
formation theories at these redshifts. 

Throughout this paper, we assume Hq = 50 km s~^ Mpc~^, and qo = 0.5, unless otherwise 
stated. At the redshift of the cluster, the luminosity distance is 5.56 hgg Gpc, the angular size 
distance is 1.70 h^Q Gpc, and the scale is 8.24 h^Q kpc per arcsec. 



2. Observations 

The discovery data of RXJ 1716.6+6708 are reported in Henry et al. (1997). Here we present 
the optical spectroscopy performed at Keck-II and the X-ray follow-up observations acquired with 
the HRI on ROSAT and with the ASCA satellite. We first describe the optical spectroscopy which 
gives a very high value for the velocity dispersion. Then we present the ASCA spectra which 
show that RXJ 1716.6+6708 has a temperature commensurate with its bolometric luminosity but 
not with its velocity dispersion. Finally, we describe in detail the HRI data and present an X-ray 
image revealing a morphology which is elongated in the same direction as the optical galaxies. 



2.1. KECK slit-mask Spectroscopy 

The region of RXJ 1716.6+6708 was observed spectroscopically on July 10 and 11, 1997, 
with the Low Resolution and Imaging Spectrograph (LRIS, pke et al. 199^ ) and slit-masks on 
the Keck-II telescope. The 300 1/mm grating blazed at 5000 A combined with our grating angle 
yielded a wavelength coverage of approximately 4500 A to 9500 A, and a pixel size of 2.45 A. A 
GG495 filter was used to eliminate the overlapping second order spectrum so as not to have any 
contamination blueward of 9500 A. Given the fact that galaxies are extended objects even at this 
redshift, we have used a slit of 1.5" in width, which gives a reduced spectral resolution of 16 A. 
For the selection of the objects to observe, we have used a very deep 3 color (B, R and I) image 
kindly supplied by Doug Clowe and published in Clowe et al. (1998). The UCSCLRIS package 
designed by Drew Phillips and collaborators at Lick Observatory was used to prepare the slit-mask 
files. We designed four masks at different position angles in order to cover the central part of the 
cluster. The four masks were exposed for 10,350s, 10,800s, 13,500s and 10,300s respectively. After 
each mask, exposure flat fields and arc calibrations were also obtained. The data were reduced 
using the standard IRAF package routines for 2-D spectra images. 

Some of the objects were observed more than once in different masks. In the end, out 



of 62 unique objects observed, 37 proved to be cluster member galaxies. The galaxies in the 
cluster are marked in Figure Q and listed in Table 1. For each galaxy the (J2000) coordinates, 
measured velocity, la error, and redshift are given. In the last column the main spectral features 
are also noted. The velocity histogram for the 37 member galaxies is shown in Figure |2[ Using 
the Timothy Beers' program, ROSTAT, we obtain an average redshift and velocity dispersion 
(corrected to first order for cosmological redshift by dividing by 1+z) which are z=0.8090±0.0051 
and aios =1522^^5q km s~^. Using different methods ( Beers, Flynn and Gebhardt 1990| , Danese 



|et al. 1980| ) we obtain the same velocity dispersion within the uncertainties. This high velocity 
dispersion was already noted by Henry et al. 1997 with only 12 galaxy redshifts (the 37 cluster 
members presented here include the 12 cluster galaxies of Henry et al. 1997). The old and new 
results are consistent with each other within the uncertainties. We will discuss the implications of 
the high velocity dispersion for this cluster in Section 3. 



2.2. ASCA Spectroscopy 



We observed RXJ 1716.6+6708 with ASCA ( p?anaka et al. 19941) on June 9 and 11, 1996 for 



115,779 and 53,191 seconds respectively. The ASCA instrumentation consists of four independent 
detectors, two Solid-State Imaging Spectrometers (SIS) sensitive in the range 0.5—9 keV (140 eV 
resolution at E=6 keV), and two Gas Imaging Spectrometers (CIS) with poorer energy resolution 
but with some efficiency up to 11 keV. The SIS observations were performed in 2-CCD Faint mode 
for High bit rate and Bright mode for Medium bit rate observations. The Faint mode data are 
converted into Bright mode data on the ground in order to combine them with the on-board Bright 
mode data. Thus all data analysis is performed on Bright data. We used the eight independent 
data-sets to extract the data in the regions of interest and then combined the datasets in 2 groups, 
one for all the CIS and one for all the SIS, for the spectral analysis. 

Data preparation was done using the software package XSELECT in FTOOLS (version 
v4.1) which allows selection of valid time intervals and removal of hot and flickering pixels. We 
used the newly reprocessed data of RXJ 1716.6+6708 as suggested by the ASCA processing 
team. The reprocessing was necessary since incorrect orbit files were used in the first processing. 
The orbit files contain a time history of ASCA's orbital parameters and the use of out-of-date 
files could cause reduced signal to noise and inaccurate fluxes. For the effective area and point 
spread function files we used the latest versions available (xrt_ea_v2_0.fits and xrt_psf_v2_0.fits). 
Clean X-ray event lists were extracted using a magnetic cut-off rigidity threshold of 6 GeV c~^ 
and the recommended minimum elevation angles and bright Earth angles to reject background 
contamination (Day et al. 1995). Events were extracted from within circles of radii 2.5 arcmin 
(SISs) and 6.25 arcmin (GISs). We rejected about 98% of cosmic ray events by using only SIS 
chip data grades 0, 2, 3 and 4. Light curves for each instrument were visually inspected to 
exclude, if present, intervals with high background or data dropouts. The individual SIS response 
matrices were generated using "sisrmg" (v 1.1, April 97 version) which takes into account the 



large time-dependence of various aspects of the SIS instruments responses. The SIS data were 
rebinned in the standard way to 512 channels applying Bright21inear with the lowest 17 channels 
(< 0.6 keV) flagged as bad. To convert SIS PHA (Pulse Height Analyzer) values to PI (Pulse 
Invariants) bins the file sisph2pi_l 10397. fits was used. The data were then regrouped so that no 
energy bin had fewer than 20 counts. Local background estimates were taken rather than the 
deep background images supplied by the ASCA Guest Observer Facility (GOF). Even if the blank 
sky background subtraction gives a better statistical quality, in the local background subtraction 
both the particle and the galactic diffuse background are properly accounted for. The background 
spectra of RXJ 1716.6+6708 were extracted with the same filters and parameters used to extract 
the source data. Background counts were extracted from within annuli centered on the source and 
with an inner radius of 6.25 (2.5) arcmin and outer radius of 11.25 (4) arcmin for the GISs (SISs). 
The details for both observations are given in Table 2. 

Spectral analysis was performed using XSPEC (vlO.O) from the software package XANADU 
(Arnaud 1996). As mentioned earlier, given the low statistics of our eight datasets (4 data-sets 



for each observation) we fitted the spectral data, together with their respective response files, 
after grouping them into 2 groups (one for the four GISs and one for the four SISs) in order to 
allow for different values of the parameters fitting the GIS and SIS data. The GIS data were 
fitted in the 0.8—9 keV band and the SIS data in the 0.6—8 keV band. The spatial resolution of 
ASCA along with the low signal-to-noise of the data permitted only a global single-temperature 



( Raymond fc Smith 1977 ) model to be fitted to the data. For the GIS dataset we fixed the 
hydrogen column density to the Galactic value along the line of sight at the cluster position, 
Nh = 3.7 X lO^'^ cm~'^ ( Dickey and Lockman 199C| ) while it was left free for the SIS dataset. 



The fitting procedure gave an anomalous measurement of an order of magnitude higher value 
for the Njif, which may reflect problems still present in the SIS calibration (see the web page 
http://heasarc.gsfc.nasa.gov/docs/asca/cal_probs.htm]! ). The normalizations for the four GIS (SIS) 



spectra were constrained to be the same. The parameters varied to provide a flt to 8 spectra were 
therefore the 2 independent normalizations, the HI value, a single emission-weighted temperature, 
and the iron abundance. A best-fit cluster rest frame temperature of kT = 5.66~!iQ^g keV was 
measured for the cluster gas and a value A= 0.43^q'2i for the metallicity. All the errors are at the 
68% confidence levels (1 a), and all the fits have an acceptable reduced x^=l-05— 1.09. Figure |^ 
shows the two-parameter x^ contours for the cluster fractional metallicity and temperature. 
The derived unabsorbed fiux in the ASCA band is f2-iofcey=(l-78±0.11)xlO^^^ erg cm~^ s^^ 
corresponding to a luminosity in the cluster rest frame of L2-ioA;eV =(8.19±0.43)xl0^^ h^p erg 
s~^, and to a bolometric luminosity Lboi=(17.40ib0.91)xl0^^ h^Q erg s~^. Converting to the 
0.5-2.0 keV ROSAT band for comparison with the ROSAT HRI data and with the ROSAT 
All-Sky survey data, and using the same assumptions as above, we derive an unabsorbed flux in 
0.5—2.0 keV equal to fo.5-2fceV =(1-66+0.09) xlO~^^ erg cm~^ s~^, and a luminosity in the cluster 
rest frame ^o.5-2keV =(4.57+0.24) x 10^^ h^Q erg s~^. This is a little higher than the fiux reported 
in Henry (1997) of (1.13+0.31) xlO^^^ erg cm~2 g-i obtained from the RASS PSPC data. 



2.3. ROSAT Imaging 

X-ray observations of RXJ 1716.6+6708 were obtained with the ROSAT High Resolution 
Imager (HRI, |Triimper 1983| ) in a single pointing in May 1997, for a net live time of 170,853s. The 
HRI operates in the ROSAT 0.1—2.4 keV energy band and provides an angular resolution of ~ 4" 
(FWHM). 

We have extracted the total counts from the cluster in a circle of 125" radius (~1.03 h^Q Mpc 
at the cluster redshift) centered at the position q(J2000) = 17^^16'" 48!2 and (5(J2000) = 67°08'22'.'4. 
This position is the centroid of the X-ray emission as determined by the maximum likelihood 
method used in the IRAF/PROS (v 2. 5. pi) task xplot.xexamine. The position of the X-ray 
centroid is fully consistent with the optical position of the cD galaxy (^^20 in Table 1) within the 
uncertainties of the HRI. The displacement in right ascension is 1.7" to the west and in declination 
is 0.4" to the north of the brightest cluster galaxy. Background counts were taken from three 
different regions free of sources within circles of 240" radius. There are 788 it 117 net counts in 
0.1—2.4 keV in the cluster region. Assuming for the temperature and the iron abundance of the 
gas the values kT=5.7 keV and A=0.4, as determined from ASCA data, and a hydrogen column 
density along the line of sight Nh=3.7x10'^^ cm~^, we determine an unabsorbed flux in 0.5—2.0 
keV of fo.5-2.ofceV =(1.42±0.21) x lO^^^ erg cm^^ g-i^ ^j^g X-ray luminosity in the 0.5-2.0 keV 
cluster rest frame is ^o.5-2keV =(3.94+0.58) x 10^'* h^Q erg s^^. Both the flux and luminosity are 
in excellent agreement within the uncertainties with the ASCA determined values (see Table 3 for 
a summary of fluxes and luminosities in the various energy bands). 

The high on-axis resolution of the HRI allows us to study in more detail the morphology of 
the X-ray emission from RXJ 1716.6+6708. To take full advantage of the intrinsic resolution of the 
instrument we created an HRI image with 2.5 x 2.5 arcsec^ pixels using software kindly provided 
by Steve Snowden. In addition to filtering the raw event list, the code also produces maps of 
the expected HRI particle background and the effective exposure time. In Figure § we show the 
X-ray iso-intensity contours from the HRI observation overlaid on an optical CCD image of the 
cluster. The HRI counts image was smoothed with a Gaussian filter of 5 arcsec width (Icr) and 
then corrected for particle background events and exposure time variations. 

As shown in Figure §, the HRI reveals that RXJ 1716.6+6708 is morphologically complex at 
X-ray wavelengths too. The emission is clearly extended over a scale of roughly 1 arcmin (~ 500 
h^Q kpc), and its shape is indicative of the non-regularity of this cluster. While the strongest 
emission originates from an almost circular peak whose location is consistent with the position 
of the cD, the lower flux contours show a pronounced centroid shift and become significantly 
elongated in the same NE-SW direction as the distribution of the galaxies^. The elongated shape 
of the emission coincides qualitatively with the dark matter contours determined by Clowe et al. 
(1998), shown in the bottom panel of Figure Q. 



^The faint source about 50 arcsec to the SE of the cluster centre is most Ukely not related to the cluster. 



2.4. X-ray Profile 

Even if /3 models applied to this cluster are poorly constrained, we will attempt in the 
following to extract an X-ray profile and derive an X-ray mass to compare to the weak lensing mass 
measured by Clowe et al. (1998). A radial profile of the X-ray emission is obtained by summing 
the HRI counts in concentric annuli of variable width centered on the peak of the X-ray emission 
and dividing by the area of the annuli. This is possible out to a radius of 125" (~1.03 h^Q Mpc at 
the cluster redshift) at which point the profile becomes indistinguishable, within the noise, from 
the background level. The profile is then fit with a /3 model ( pavaliere and Fusco-Femiano 1976| ) 
described by the standard form S{r) = So{l + r^/r^)"^^"*"-^/^, with So = central surface brightness, 
Tc = core radius and f3 = slope parameter. The values of the best fit (90% confidence level) are 
So = 3.24 X 10^13 erg cm"^ g-i arcmin"2, r^ = 6'.'8 ± 6'.'1 and (3 = 0.42+[];o9. Both the core radius 
and the /3 parameter are poorly constrained given the irregular X-ray morphology. At the distance 
of the cluster this core radius corresponds to 56 h^Q kpc with a similar size error (from 5.7 to 
106 kpc). The size of the core radius could indicate the presence of a cooling flow. However, 
there is no evidence for excess emission in the very core of the cluster, thus any cooling flow must 
be of marginal significance. The existence of a cooling flow in RXJ 1716.6+6708 would also be 
inconsistent with the emerging view that cooling flows form preferably in dynamically undisturbed 
clusters (see Peres et al. among others). A fit of a point source plus the f3 model to the data 
results in a contribution of less than 15% of the point-like component. To confirm or disprove the 
presence of a central excess a spectral-spatial fit would be required. This is not possible with the 
X-ray data in hand, we have nonetheless attempted to fit a 2-temperature model to the ASCA 
spectra. The addition of the second temperature model does not improve significantly the fit with 
respect to the 1-Temperature model (F-test probability less than 2%). 

Small core radii (less than 40-50 kpc) and small f3 values (~ 0.5), such as this cluster seems 
to possess, have been observed in other clusters; see for instance MS1512+3647 {re = 6'.'9 it I'.'l 
(42.6±6.7 h^o^kpc), /3 = 0.524 ± 0.031; [Hamana et al. 1997|) or MS0440+0204 (r^ = 6'.'4 ± I'.'l 



(26.7±4.5 h^o^ kpc), (3 = 0.45 ± 0.03; |Gioia et al. 1998| ). The small size of the core radfi and /3s 



for these X-ray selected clusters is an issue that will be investigated in future work. Here we use 
these two correlated parameters to derive an estimate of the cluster mass. 



2.5. Mass estimate: X-ray 

The mass determined using X-ray data does depend on assumptions involving spherical 
symmetry and hydrostatic equilibrium, thus some of the assumptions may not apply in this case. 
These assumptions are found to be valid on the average in N-body simulation studies of cluster 
formation by Schindler (1996) and Evrard et al. (1996). Moreover, from extensive simulations 
Evrard et al. (1996) have shown that errors on cluster masses, using scaling relations between 
mass and temperature, are typically of the order of 20% even in cases, like RXJ 1716.6+6708, 



where some substructure is present. Allen et al. 1996 caution against the use of single-phase 



analyses with ROSAT data, since the presence of a cooling flow with distributed mass deposition 
implies that the central intracluster medium (ICM) has a range of temperatures and densities at 
any particular radius: i.e. the ICM is multi-phase. There is no evidence of any significant cooling 
flow in RXJ 1716.6+6708, thus we feel that the assumptions adopted for the deprojection outlined 
below are largely correct. We should bear in mind though that the derived quantities using the 
HRI profile should be interpreted with some caution given the asymmetry of the X-ray contours 
map that suggest that RXJ 1716.6+6708 is not a well relaxed cluster. 

From the HRI surface brightness profile, assuming a constant temperature, we have derived 
the three-dimensional density distribution of the gas from the two-dimensional image. The 
integrated mass in all forms can then be derived as a function of radius directly from the equation 
of the hydrostatic equilibrium 
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The radial dependence of the gas density pg is given by the ROSAT HRI observations. Tg is 
the intracluster temperature, // is the mean molecular weight of the gas, and rup is the proton 
mass. A constant intracluster gas temperature of 5.66 keV as measured by ASCA is assumed. 
Our estimates for the gas and gravitational mass within 125" (~1.03 h^Q Mpc at the cluster 
redshift) are (8.9+2.1) x lO^^ Yu^^^'^Mq and (2.8+0.3) x 10^"^ h^o^M© . The corresponding gas 
mass fraction, Mgas/Mtot=(31+8)% h^Q , is higher than found by other investigators (11—25% 







hgQ David et al. 1995] and [White and Fabian 1995). More recently Evrard (1997) has published 



a revised average value of 17% for the gas fraction using the above referenced compilations. There 
is not yet a consensus on the value of the hot gas to total mass ratio in clusters, and on the 
fact whether this ratio is really representative of the fraction of baryons in the universe. The 
measurements seem to indicate that the gas fraction may vary from cluster to cluster and may 
even vary in the same cluster as a function of radius (see Donahue 1998). 

We compare now the total mass from X-rays to the mass derived from the gravitational shear 
signal by Clowe et al. (1998). These authors find in roughly the same region (their 500 h~^kpc 
aperture corresponds to 1 Mpc in our cosmology) a value of (5.2+1.8) x 10^^ h^Q Mq . The weak 
lensing mass is computed assuming the background galaxies lie on a sheet at z=2 and is higher, 
even within the uncertainties, than the X-ray derived mass. We note however that the Clowe et 
al. (1998) mass determination takes into account also the second clump of matter associated with 
the NE group of galaxies which is clearly separated in their mass distribution map from the main 
cluster mass (see Figure 0). The X-ray gravitational mass comes only from the main optical clump 
centered on the cD. In addition the weak lensing center of mass of the main cluster is displaced 
with respect to the center of the light and X-ray peaks of the cluster; it is located 27" east and 
11" north of the cD (see Figure |^ but consistent within the uncertainties of the weak lensing 
positions. The weak lensing mass distribution also resembles a filament of structure elongated in 
the same direction as the chain of optical galaxies. 



3. Results and Discussion 

In the following discussion we will focus on the implications of measured quantities such 
as the high velocity dispersion, or low temperature, rather than derived quantities such as the 
gravitational mass from the X-ray data. First we would like to note that using three completely 
different instruments (ROSAT PSPC, ASCA and ROSAT HRI) we obtain within the uncertainties 
a consistent flux for the cluster. This is more remarkable than it seems since one has to keep in 
mind that the discovery data of RXJ 1716.6+6708 consisted of only 33 net counts. 

There are not many X-ray selected clusters known at this high redshift; the other published 
examples are the 2 clusters from the EMSS ( |Gioia et al. 1990D , MS1054-03 at 0.83 and 
MS1137+66 at 0.78 ( pioia and Luppino 1994 ). Donahue et al. (1998) have made a follow-up 



study of MS1054— 03 similar to the one presented here. We will compare in what follows the 
similarities and differences between RXJ 1716.6+6708 and MS1054— 03. A more detailed study of 
the other Medium Survey cluster, MS1137+66, has been submitted (Donahue et al. 1999). 

RXJ 1716.6+6708 and MS1054-03, both at a redshift higher than 0.8, show a filamentary 
appearance in their optical morphology. The Medium Survey cluster has an optical elongation 
in the east-west direction with the HRI contours extended over a region of about 1' in the same 
direction as the optical galaxies ( ponahue et al. 1998|) . RXJ 1716.6+6708 also has an extended 
X-ray emission of about 1', with the center of the X-rays on the cD galaxy and the elongation in 
the same direction as the optical galaxies. 

Similarly to MS1054— 03, RXJ 1716.6+6708 has a high velocity dispersion which can be 
interpreted as a signature of non-virialization. However, differently from MS1054— 03, we are 
dealing here with a rather cool cluster (kT=5.66 keV against kT=12.3 keV). The temperature of 
RXJ 1716.6+6708, 5.66^g'5g, is commensurate with the predictions from its X-ray luminosity from 
the Lx— Tx relation by David et al. (1993) and Arnaud and Evrard (1998). There is an extensive 
literature on the correlation between these two basic and measurable quantities (Edge & Stewart 
1991, Ebeling 1993, David et al. 1993, Fabian et al. 1994, Mushotzky & Scharf 1997, Markevitch 
1998, Arnaud and Evrard, 1998). Comparing the bolometric luminosity of RXJ 1716.6+6708 with 



the best fit relation, log(Lx)=(2.88±0.15) log(T/6keV) +(45.06+0.03) obtained by |Arnaud and 



Evrard 1998| , analyzing a sample of 24 low-z clusters with accurate temperature measurements 



and absence of strong cooling flows, we would expect for RXJ 1716.6+6708 a temperature of 6.9 
keV. We measure 5.66lQ'5g. The Lx—T^x relation does not seem to evolve much with redshift since 
z=0.4 (Mushotzky & Scharf 1997; Henry 1997) and possibly since z=0.8 (see study on MS1054-03 



by Donahue et al. 1998). Modest evolution is predicted in some evolution models with JIq =1 
if the gas is pre-heated before falling into the cluster potential (Evrard and Henry, 1991; Kaiser, 
1991). 

Next we check the a—Tx relationship. The measured velocity dispersion of RXJ 1716.6+6708, 
1522]'^^5Q km s^^, is much higher than expected from its temperature. A large number of authors 
(see Table 5 in Girardi et al., 1996, or Table 2 in Wu, Fang and Xu, 1998, for an exhaustive list 



10 



of papers on the subject) have attempted to determine the a—Tx using different cluster samples 
in order to test the dynamical properties of clusters. Using the derived value from the HRI data 
of /3=0.42 and the measured temperature of 5.7 keV from the ASCA data for the temperature 

2 

of the gas in the equation d = ^'^^'^" we obtain a velocity dispersion of about 600 km s~^. Even 

"'-' gas 

assuming energy equipartition between the galaxies and the gas in the cluster (/3=1) we still 
obtain a velocity dispersion much lower than measured (~ 929l4g km s~^ vs 1522 km s^^). In 
other words in an isothermal potential the mass distribution of the cluster does not correspond 
to the isotropic one-dimensional velocity dispersion. This is rather at odds with results of hot, 
high-z clusters recently observed. The velocity dispersions of three high (>0.5) redshift clusters, 
MS1054-03, MS0016+16 and MS0451-03 (see discussion in Donahue et al. 1998) ah reflect their 
temperatures according to the above equation with /3=1. Girardi et al. (1996) have derived a 
best fit relation between the velocity dispersion and the X-ray temperature, with more than 30% 
reduced scatter with respect to previous work (Edge and Stewart 1991; Lubin and Bahcall 1993; 
Bird, Mushotzky and Metzler, 1995; Wu, Fang and Xu 1998, among others). Girardi et al. (1996) 
have taken into account distortions in the velocity fields, asphericity of the cluster or presence of 
substructures to derive their best fit relation, log(<T)=(2.53±0.04)-|-(0.61±0.05)log(T). If we plug 
in the temperature of RXJ 1716.6-1-6708 in the above relation, the resulting velocity dispersion 
value is discrepant with what is measured. A temperature of 11.5^^g keV would be expected for 
<7=1522 kms~^ 

RXJ 1716.6-1-6708 may be an example of cluster which has not reached virial equilibrium, its 
dynamical state may be in large part dominated by infall or merging and consequently the velocity 
dispersion is not representative of the virial temperature of the cluster. If RXJ 1716.6-1-6708 is 
actually composed of two or more distinct gravitational components the situation should be 
apparent from the velocity histogram (which is Gaussian from Fig. 2), unless the separation of 
the cluster components is so small to be statistically of little significance. The high value of the 
velocity dispersion can be due to a number of causes. If there is any fraction of infalling galaxies 
which are bound to the cluster but not yet virialized, they could inflate the velocity dispersion. 
These galaxies in the cluster would be moving on radial orbits. To see if any substructure is 
present in the redshift space we have considered the galaxies in Figure || lying outside a circle (0.5 
Mpc radius) centered on the cD, to the north-east and to the south-west of a dividing line passing 
through the cD with position angle of about -|-135° East of North. There are eleven galaxies to the 
NE (labeled in Table 1 as 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 and 25) and nine galaxies to the SW (labeled 
27, 28, 29, 30, 31, 32, 33, 34 and 37). The mean redshift of the galaxies to the NE is 0.8107±0.0018 
(1 a on the mean), and the mean redshift of the galaxies to the SW is 0.8021ib0.0027, which 
implies a difference in redshift between the two components of Az=0.0086±0.0032 (1425 km s~^ 
in velocity space), a small effect that is significant at 2.65 a. It might well be that galaxies are 
infalling towards the center of the cluster, and while some galaxies have already reached (or 
crossed) the core region, some others are still moving along radial orbits towards its center. 

We recall here that numerical simulations have shown that the universe is composed of a 
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web of filaments and voids on the scale of hundreds of niegaparsecs. The clusters of galaxies 
form at the intersections of these filaments and grow with time from quasi-spherical systems to 
the spherical objects that we observe today as gas, galaxies and dark matter potential wells. 
Filamentary structures have been observed in less distant clusters (i.e. A2125, Wang et al. 1998; 
the ABCG 85/87/89 complex, Durret et al. 1998). The initial formation of protoclusters is often 
described as matter flowing along filaments ( [Bond, Kofman and Pogosyan, 1996| ) with the X-ray 
flux maximum at the impact point of two colliding streams of matter. In this scenario the hot 
gas and the galaxies would not be in hydrostatic equilibrium and the transient velocity dispersion 
could be higher than expected from a virialized system. We are probably witnessing this process 
in this very distant X-ray selected cluster, and RXJ 1716.6-1-6708 might be more properly called a 
protocluster. 
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to Timothy Beers for the use of his ROSTAT program to calculate velocity dispersions. Steve 
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AST95-00515, from NASA-STScI grants GO-5402.01-93A and GO-05987.02-94A, from NASA 
grants NAG5-1880 and NAG5-2523, from SAO contract SV4-64008 and from CNR-ASI grants is 
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Fig. 1. — The image is a 1024x1024 subarray extracted from the center of a 4500s exposure in the 
I-band taken by Luppino and Metzger with the UH 8Kx8K CCD mosaic-camera on the CFHT 
prime focus. The exposure was made up of 5 dithered pointings each 900s long. This image has a 
seeing of 0" .7 FWHM and spans 7'. 2x7'. 2 (3.5x3.5 h^Q^Mpc at the redshift of the cluster). North 
is up and East to the left. 



Fig. 2. — Distribution of velocities for the entire sample of 37 galaxies. The curve shows a Gaussian 
distribution with width given by the derived aios and normalized to the 37 galaxies considered 
cluster members. The K-S test plus a suite of other robust diagnostics indicate the distribution is 
Gaussian with a probability >90%. 



Fig. 3. — Two-dimensional x^ contours at 68.3% and 90% confidence levels {Ax^ =2.30 and 4.61) 
for the cluster iron abundance in units of solar abundance and temperature in keV. These contours 
represent the confidence contours of simultaneous fits to GIS and SIS spectra binned such that each 
bin contains a minimum of 20 counts. 



Fig. 4. — The upper panel shows a CFHT I-band image of the central region of 
RXJ 1716.6-1-6708 with contours of the HRI flux overlaid. The HRI raw counts were smoothed 
with a Gaussian of 5 arcsec width (la) and then corrected for particle background contamination 
and exposure time variations. The remaining X-ray background in this image is 4.3 x 10^^ ct s^^ 
arcsec"^; the contours are placed at 1.8 x background x 1.15". The bottom panel shows a Keck-II 
R-band image of the same area with the weak lensing contours from Clowe et al. (1998) overlaid. 
Either panel covers an area of 2.75 x 2.75 arcmin^ (1.3 x 1.3 h^Q Mpc^ at the redshift of the cluster). 



Table 1. Spectroscopic data for the member galaxies of RXJ 1716.6+6708 



ID 


a(J2000) 


5(J2000) 


cz (km/s) 


cz (err) 


z 


Spectral Features 


1 


17:17:19.5 


+67:11:48.9 


241309 


150 


0.8049 


H+K, CN 


2 


17:17:14.7 


+67:11:35.3 


244127 


900 


0.8143 


[Oil] 


3 


17:17:14.2 


+67:11:39.7 


241699 


60 


0.8062 


H+K, CN 


4 


17:17:09.7 


+67:10:48.9 


244607 


240 


0.8159 


[Oil] 


5 


17:16:59.1 


+67:10:18.6 


246046 


150 


0.8207 


H+K, CN, H7 


6 


17:17:08.9 


+67:09:27.1 


243857 


270 


0.8143 


H+K 


7 


17:17:01.0 


+67:09:32.1 


244187 


120 


0.8145 


H+K, G-band 


8 


17:17:01.2 


+67:09:24.0 


242418 


90 


0.8086 


H+K, CN 


9 


17:17:03.8 


+67:08:59.6 


242118 


90 


0.8076 


[Oil], H+K, m 


10 


17:16:56.9 


+67:09:04.6 


240110 


90 


0.8009 


H+K, CN 


11 


17:16:55.1 


+67:09:05.1 


240350 


90 


0.8017 


H+K, CN, H(5 


12 


17:16:51.0 


+67:08:51.4 


246016 


90 


0.8206 


H+K, CN, G-band 


13 


17:16:54.5 


+67:08:42.5 


239630 


30 


0.7993 


H+K, CN 


14 


17:16:55.0 


+67:08:32.7 


246315 


900 


0.8216 


[Oil] 


15 


17:16:51.4 


+67:08:27.6 


241639 


160 


0.8060 


H+K, CN 


16 


17:16:54.1 


+67:08:21.9 


244817 


900 


0.8166 


Call-break 


17 


17:16:53.6 


+67:08:07.6 


240199 


900 


0.8012 


H+K 


18 


17:16:50.1 


+67:08:20.7 


238191 


180 


0.7945 


H+K, CN 


19 


17:16:50.1 


+67:08:22.4 


239240 


900 


0.7980 


[Oil] 


20 


17:16:48.5 


+67:08:22.0 


247515 


450 


0.8256 


cD GALAXY, H+K 


21 


17:16:47.3 


+67:08:17.6 


243467 


180 


0.8121 


H+K, CN, H5 


22 


17:16:49.1 


+67:08:22.0 


246106 


240 


0.8209 


H+K, CN 


23 


17:16:46.7 


+67:08:16.0 


243078 


120 


0.8108 


H+K, CN 


24 


17:16:44.4 


+67:08:11.0 


247095 


150 


0.8242 


Call-break, CN 


25 


17:16:43.7 


+67:09:43.9 


242838 


900 


0.8100 


[Oil] 


26 


17:16:39.0 


+67:08:33.0 


242868 


30 


0.8101 


[Oil], H+K, US 


27 


17:16:36.9 


+67:08:28.9 


238251 


60 


0.7947 


[Oil], H+K, CN, G-band, 4C 67.26 


28 


17:16:36.9 


+67:08:14.1 


241039 


90 


0.8040 


[Oil], CN, H+K, H7 


29 


17:16:34.4 


+67:08:02.7 


242118 


900 


0.8076 


[Oil] 


30 


17:16:23.5 


+67:07:43.4 


242598 


900 


0.8092 


[Oil] 


31 


17:16:32.6 


+67:07:36.4 


244157 


150 


0.8144 


[Oil], H+K, G-band 


32 


17:16:35.7 


+67:07:48.2 


239630 


90 


0.7993 


[Oil], CN, H+K 


33 


17:16:37.6 


+67:07:30.4 


241159 


120 


0.8044 


H+K 


34 


17:16:42.6 


+67:06:42.9 


237771 


900 


0.7931 


[Oil] 


35 


17:16:49.0 


+67:07:20.3 


242928 


60 


0.8103 


H+K, CN, G-band, H7 


36 


17:16:54.7 


+67:07:38.5 


246945 


300 


0.8237 


H+K, CN, G-band 


37 


17:16:58.5 


+67:06:52.2 


237561 


30 


0.7924 


H+K 



Table 2. ASCA observing information 



Instrument 


Arcmin Radius 


GIS2-obsl 


6.25 


GIS2-obs2 


6.25 


GIS3-obsl 


6.25 


GIS3-obs2 


6.25 


SISO-obsl 


2.5 


SIS0-obs2 


2.5 


SISl-obsl 


2.5 


SISl-obs2 


2.5 



Cts s ^ Screened Exposure in s 



(3.618±0.750)xl0"^ 54,713 

(5.471±1.154)xl0"^ 25,109 

(7.604±0.824)xl0"^ 54,713 

(8.781±1.214)xl0"^ 25,105 

(5.171±0.596)xl0"^ 45,687 

(6.278±0.865)xl0"^ 21,005 

(4.691±0.619)xl0"^ 45,694 

(5.393±0.905)xl0-^ 20,981 



Table 3. X-ray data quantities 



Instrument 


f 0.5-2 keV 


f2-io keV 


Lo.5-2 keV 




L2-10 keV 


LsoL keV 




erg cm^^ s^^ 


erg cm^^ s^^ 


erg s~^ 




erg s"^ 


erg s"^ 


PSPC 


(1.13±0.31)10-^^ 


_ 


(3.20±0.90)10*'' 




_ 


_ 


ASCA 


(1.66±0.10)10-i^ 


(1.78±0. 11)10"^^ 


(4.57±0.24)10'''' 


(8 


.19±0.43)10*'' 


(17.4±0.91)10*^ 


HRI 


(1.42±0.21)10-^^ 


- 


(3.94±0.58)10^'' 




- 


- 



This figure "Gioia.figl.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/9902277vl 



This figure "Gioia.fig2.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/9902277vl 



This figure "Gioia.fig3.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/9902277vl 



This figure "Gioia.fig4.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/9902277vl 



